Search for pair-production of strongly-interacting particles 
decaying to pairs of jets in pp collisions at a/s = 1.96 TeV 
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We present a search for the pair-production of a non-standard-model strongly-interacting particle 
that decays to a pair of quarks or gluons, leading to a final state with four hadronic jets. We 
consider both non-resonant production via an intermediate gluon as well as resonant production via 
a distinct non-standard-model intermediate strongly-interacting particle. We use data collected by 
the CDF experiment in proton-antiproton collisions at y^i — f .96 TeV corresponding to an integrated 
luminosity of 6.6 fb~^ . We find the data to be consistent with standard model predictions. We report 
limits on cr(pp — > jjjj) as a function of the masses of the hypothetical intermediate particles. Upper 
limits on the production cross sections for non-standard-model particles in several resonant and 
non-resonant processes are also derived. 
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One of the few hints of possible physics beyond the 
standard model (SM) at the TeV scale is the anomalous 
top-quark forward-backward asymmetry Afb observed at 
the Tevatron [TH3]. This asymmetry could be generated 
by non-SM physics through the production of top-quark 
pairs via a light axi-gluon H], a particle with axial cou- 
plings to quarks, that interferes with standard model 
(SM) tt production to produce the observed asymme- 
try. The axi-gluon would be visible in its alternate decay 
mode to low-mass strongly-interacting particles, each of 
which decays to a pair of jets [5] yielding a four-jet final 
state. This final state is of broad interest, as various mod- 
els predict pair-production of strongly-interacting par- 
ticles decaying to jet pairs with no intermediate reso- 
nance [HI E] and i?-parity- violating supersymmetric the- 
ories [5] predict pair-production of light partners of the 
top quark (stop quarks), each decaying into to pairs of 
light quarks. 

The masses of the axi-gluon and its strongly- 
interacting decay products are not predicted, but must be 
fairly light (< 400 GeV/c^) to explain the Afi, measure- 
ment . The LHC experiments have excellent sensitivity 
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FIG. 1: Diagrams for resonant (left, via X) and non-resonant 
(right) pair-production of Y particles, with subsequent decays 
to pairs of gluons. Other models, with final-state quarks, are 
also considered. 



at high mass due to the large center-of-mass energy, but 
difhculties at low mass due to high background rates. 
The ATLAS experiment ruled out masses between 100 
and 150 GeV/c^ [ID]; CMS ruled out masses between 250 
and 740 GeV/c^ [TT]. No experimental bounds exist for 
such non-SM particles with masses below 100 GeV/c^ for 
non-resonant pair-production of di-jet resonances; there 
are no current limits on resonant production. 

In this Letter we report a search for both non-resonant 
and resonant production of pairs of strongly-interacting 
particles, each of which decays to a pair of jets. Rather 
than probing a specific theory, we construct a simplified 
model with the minimal particle content. In the non- 
resonant case, we consider the production process pp — >■ 
YY -^ii jj, with the mass of the hypothetical Y state, 
my as a single free parameter. In the resonant case, 
— ^ X — > YY -^jj jj, we also explore the mass of the 
X state, mx (Fig. [T]). In both cases, we assume that the 
natural width of the particles is small compared to the 
experimental resolution. 

We analyze a sample of events corresponding to an in- 
tegrated luminosity of 6.6±0.5 fb~^ recorded by the CDF 
II detector ^2 , a general purpose detector designed to 
study pp collisions at -/s = 1.96 TeV produced by the 
Fermilab Tevatron collider. The tracking system consists 
of a silicon microstrip tracker and a drift chamber im- 
mersed in a 1.4 T axial magnetic field [T3]. Electromag- 
netic and hadronic calorimeters surrounding the tracking 
system measure particle energies, with muon detection 
provided by an additional system of drift chambers lo- 
cated outside the calorimeters. 

We reconstruct jets in the calorimeter using the jet- 
CLU [13] algorithm with a clustering radius of 0.4 in 77 — 
space [15j . and calibrated using the techniques outlined 
in Ref. [IH] . Events are selected online (triggered) by the 
requirement of three jets, each with Et > 20 GeV and 
with Ejots-E'T > 130 GeV [15]. The data set used in this 
search is limited to 6.6 fb"""^ because the trigger selection 
was not available in early data. After trigger selection, 
events are retained if at least four jets are found with 
Et > 15 GeV and \t]\ < 2.4. 

We model resonant and non-resonant production with 
MADGRAPH5 171 version 1.4.8.4 and the CTEQ6l1 1181 
parton distribution functions (PDF). Additional parton 
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Resonance mass or [GeV/c^] 

FIG. 2: Overall efficiency, including trigger and selection 
requirements. Efficiency is shown for several simulated non- 
resonant YY -^jj jj samples with varying my- The shaded 
band shows the uncertainty. In addition, efficiency is shown 
for several simulated resonant X — >■ YY -^jj jj samples with 
varying mx and my. The uncertainty is not shown but is 
similar to the non-resonant case. The turn-on curve is deter- 
mined largely by the trigger requirement that Ejets-ET > 130 
GeV. 

radiation, hadronization, and underlying-event modeling 
are described by PYTHIA [TO] version 6.420. The detec- 
tor response for all simulated samples is modeled by the 
GEANT-based CDF II detector simulation [20] • 

The trigger and selection requirements have an effi- 
ciency on the signal up to 90% if Tj-^atsEr exceeds signif- 
icantly the 130 GeV trigger threshold. For events with 
smaller Ejets^^T, the efficiency decreases rapidly (Fig. [2|. 
In the non-resonant-production model, the Ejets^/T is 
strongly correlated with my. In the resonant-production 
model it is correlated with mx ; additionally if mx — 2my 
is large, the pr of the resulting Y is large, which leads 
to a small opening angle of its decay products and a loss 
of efficiency due to merged jets. The trigger efficiency is 
measured in simulated events, and uncertainties derived 
from validation in disjoint samples; the measured trig- 
ger efficiency and uncertainty are applied to the signal 
model. 

To reconstruct the di-jet resonance, we consider the 
four leading jets and evaluate the invariant mass of each 
of the di-jet pairs in the three permutations, choosing 
the permutation with the smallest mass difference be- 
tween the pairs. As the pair masses are correlated, we 
take the mean of the two pair masses as the estimate of 
the di-jet resonance mass. To reduce backgrounds, we 
require that the relative mass difference between the two 
pairs is less than 50%, and that the production angle 
6* of the di-jet resonance in the YY pair center-of-mass 
frame satisfies cos(6'*) < 0.9. In the resonant production 
analysis, we calculate the four-jet invariant mass. No spe- 
cific my-dependent selections are made; the requirement 
that the relative di-jet mass difference be small ensures 
compatibility with the X — YY hypothesis. Figures [3] 
and [4] show the observed di-jet and four-jet spectra, re- 
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Data PP^YY- 
m^=70 GeV/c^ o=200 pb 
mY=120 GeV/c^o=50 pb 
m^=400 GeV/c^ o=0.5 pb 
Fit 
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FIG. 3: Reconstructed mean di-jet mass in events with four 
jets. Parametric fit and several signal hypotheses overlaid in 
(a). Relative difference between the observed data and the fit 
in each bin shown in (b). 



FIG. 4: Reconstructed four-jet mass in events with four jets. 
Parametric fit and several signal hypotheses overlaid in (a). 
Relative difference between the observed data and the fit in 
each bin shown in (b). 



spectively. 

The dominant background originates from standard 
QCD multi-jet production. We model this background 
contribution using a parametric function which is fit to 
the reconstructed mass spectrum of the observed data. 
The function is a piece-wise combination of a third-order 
polynomial to describe the turn-on region, a third-order 
polynomial to describe the peak region, and a double ex- 
ponential of the form /(m) = aie~^™~°^-' ^/"^ to describe 
the falling spectrum. The parametric functional form 
was chosen to be flexible enough to describe the multi- 
jet mass spectrum, but rigid enough to avoid accurately 
describing a spectrum which includes a narrow resonance, 
so that in the presence of a narrow feature a signal-plus- 
background hypothesis would be preferred. For the di- 
jet mass, the ranges used are [35, 82.5], [82.5, 140], and 
[140,700] GeV/c^; for the four-jet mass, the ranges used 
are [115, 185], [185, 330], and [330,800] GeV/c^. The 
functional form and ranges were chosen based on their 
ability to accurately describe the mass spectra of simu- 
lated multi-jet events generated by ALPGEN [21] version 
2.10. 

The dominant source of systematic uncertainty is due 
to the multi-jet background model. The functional form 
is an approximation, which even in the absence of a nar- 
row feature may deviate from the observed spectrum. 
We estimate the impact of these potential deviations by 
measuring their magnitude in two background-enriched 
control samples. These two control samples are adjacent 
to the signal region and capture the expected deviations 
in two independent directions. The first requires a large 
relative di-jet mass difference, greater than 50%, and the 
second requires cos{9*) > 0.9. The observed relative de- 
viations are then applied to the observed spectrum in 
the signal region to estimate the magnitude of spurious 



deviations due to possible mismodeling. In addition, we 
verify that the fitting procedure gives an unbiased esti- 
mate of the signal rate. 

An additional uncertainty is due to knowledge of the 
trigger efficiency extracted from the simulated sig- 
nal samples, varying from 20% relative at Ejets-E-T = 120 
GeV to 10% above ^etsEx = 200 GeV. Uncertainties 
in the levels of parton radiation |23| and in the calibra- 
tion of the jet energy and resolution modeling [TB] also 
contribute to uncertainties in the trigger and selection 
efficiency and reconstructed mass spectrum of the signal 
samples. These uncertainties are small (< 10%) relative 
to the fitting and trigger uncertainties. 

In the non-resonant analysis, for each Y mass hypoth- 
esis, we fit the most likely value of the Y pair-production 
cross section (cff) by performing a maximum likelihood 
fit of the binned di-jet mass distribution, allowing for sys- 
tematic and statistical fluctuations via template morph- 
ing [21]. The likelihood takes the form of 

L{aYY) = Y[ /bg(o^) + '^YY C e /sig, 

bin i 

where /bg('3) is the parametric function with nuisance 
parameters a defined above to describe the background 
spectrum, fsig is a normalized template of the expected 
shape of the signal determined from simulated events, 
and Ce is the product of the integrated luminosity and 
efficiency. No evidence is found for the presence of pair- 
production of di-jet resonances and upper limits on Y 
pair-production at 95% confidence level (C.L.) are set. 

Limits are calculated using the CLs [25] method by 
repeating the measurement on sets of simulated exper- 
iments that include signal contributions corresponding 
to various hypothetical production cross-sections, and 
variation of systematic uncertainties. The values of nui- 
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FIG. 5: Upper limit on signal production rate at 95% C.L. 
Expected and observed upper limits on (7{pp — >■ YY -^jjjj) 
versus my in the non-resonant analysis shown in (a). Two 
signal hypotheses are shown, see text for details. Observed 
limits on (j{pp X ^ YY -^jjjj) versus mx and my shown 
in (b). Circles indicate the true values of the parameters used 
in each ensemble of simulated samples used to evaluate the 
limits; intermediate values are interpolated. 



sance parameters are not fit in the experiments. The 
observed limits are consistent with expectation for the 
background-only hypothesis. The resonant analysis is 
very similar, but is done as a function of the X mass 
hypothesis, fitting the four-jet mass distribution for the 
most likely value of X production cross section, <tx- 

In the non-resonant case, this analysis sets limits on 
coioron or stop-quark pair production, excluding 50-100 
GeV/c^ and 50-125 /c^, respectively; see Table [l] and the 
top of Fig. [5] The uncertainty on the theoretical cross- 
section prediction comes from two sources summed in 
quadrature. The first uncertainty is the envelope of the 
PDF uncertainties from the CTEQ uncertainties and an 
alternative PDF choice, mstw2008lo [26] (5% relative). 
The second uncertainty comes from a variation of the 
renormalization and factorization scales by a factor of 
two in each direction from their default values of the 
per-event mass scale. These theoretical uncertainties are 
illustrated in Figure [Sj 

In the resonant case, this analysis excludes axi-gluon 
(A) production leading to pairs of a particles and four- 



TABLE I: Observed and expected 95% C.L. upper limits on 
a{pp — >■ YY -^jj jj) for several values of my. Also shown 
are theoretical predictions for coioron pair production [HI [7] 
or stop-quark pair production with i?-parity-violating decay 



Mass 


Expected Observed 


Coioron 


Stop quarks 


(GeV/c^) 


(Pb) 


(pb) 


(pb) 


(pb) 


50 


240 


250 


320 


570 


70 


75 


62 


180 


100 


90 


8.2 


5.9 


62 


26 


100 


11 


17 


37 


15 


125 


14 


11 


11 


4.4 


150 


37 


46 


3.7 


1.5 


200 


4.5 


2.0 


0.60 


0.25 


250 


2.7 


1.5 


0.11 


5.4 X 10"^ 


300 


2.0 


3.0 


2.9 X 10"^ 


1.3 X 10"^ 


400 


1.1 


1.5 


1.7 X 10"^ 


7.2 X 10"* 


500 


0.3 


0.3 


8.5 X 19"^ 


3.6 X 10"^ 



TABLE II: Observed and expected 95% C.L. upper limits on 
a{pp — >■ X — > YY -^jj jj) for several values of my and mx- 
Also shown are theoretical predictions for axi-gluon produc- 



assuming 


coupling 


to quarks of Cq 


= 0.4 5, 9, 




mx 


my 


Expected Observed Axi- 


gluon 


(GeV/c^ 


) (GeV/c^ 


) (pb) 


(pb) 


(pb) 


150 


50 


641.2 


431.1 


5600 




70 


209.6 


270.6 




175 


50 


66.8 


78.9 


3500 




70 


111.5 


163.9 




200 


50 


13.8 


9.5 


2200 




70 


30.4 


91.5 






90 


17.8 


100.4 




225 


50 


18.0 


26.0 


1750 




70 


20.7 


25.0 






90 


20.9 


25.3 




250 


50 


6.2 


2.0 


1000 




70 


4.0 


3.6 






90 


5.1 


2.8 




275 


50 


6.5 


1.2 


850 




70 


7.7 


1.3 






90 


9.7 


1.4 




300 


50 


5.0 


7.1 


540 




70 


2.4 


2.6 






90 


1.7 


1.0 






140 


1.8 


1.2 




400 


50 


15.5 


6.8 


170 




70 


15.0 


20.2 






90 


30.6 


52.8 






140 


41.0 


74.6 






180 


46.9 


79.1 




500 


50 


20.7 


6.8 


60 




70 


15.9 


4.7 






90 


17.7 


5.9 






140 


25.2 


7.0 






180 


26.7 


8.0 






220 


29.7 


9.3 





7 



gluon final state for niA € [150, 400], rrio- G [50,m^/2] in 
the case of coupling to quarks Cq = 0.4 (see Table [Tl| and 
the bottom of Fig. [5| which is close to the value required 
to explain the top-quark Afi, result [S]. To be consistent 
with this analysis, the couplings would have to be smaller 
by an order of magnitude. Maintaining consistency with 
the top-quark Af}, result would require different couplings 
to light quarks and heavy quarks, with the heavy-quark 
coupling approaching the perturbative limit, Cq < 1. 

We thank Martin Schmaltz, Gustavo Tavares, Can 
Kilic, Bogdan Dobrescu, Dirk Zerwas and Felix Yu for 
useful suggestions and technical advice. We thank the 
Fermilab staff and the technical staffs of the partici- 
pating institutions for their vital contributions. This 
work was supported by the U.S. Department of Energy 
and National Science Foundation; the Italian Istituto 
Nazionale di Fisica Nucleare; the Ministry of Education, 
Culture, Sports, Science and Technology of Japan; the 
Natural Sciences and Engineering Research Council of 
Canada; the National Science Council of the Republic 
of China; the Swiss National Science Foundation; the 
A. P. Sloan Foundation; the Bundesministerium fiir Bil- 
dung und Forschung, Germany; the Korean World Class 
University Program, the National Research Foundation 
of Korea; the Science and Technology Facilities Council 
and the Royal Society, UK; the Russian Foundation for 
Basic Research; the Ministerio de Ciencia e Innovacion, 
and Programa Consolider-Ingenio 2010, Spain; the Slo- 
vak R&D Agency; the Academy of Finland; and the Aus- 
tralian Research Council (ARC). 



[1] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 

101, 202001 (2008). 
[2] V. M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 

100, 142002 (2008). 
[3] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 

83, 112003 (2011). 
[4] P. H. Frampton and S. L. Glashow, Phys. Lett. B 190, 

157 (1987). 

[5] G. Marques Tavares and M. Schmaltz, Phys. Rev. D 84, 
054008 (2011). 



[6] Y. Bai and B. A. Dobrescu, J. High Energy Phys. 07 
(2011) 100. 

[7] C. Kilic, T. Okui, and R. Sundrum, J. High Energy Phys. 
07 (2008) 038; C. Kilic, S. Schumann, and M. Son, J. 
High Energy Phys. 04 (2009) 128. 
[8] S. P. Martin, In Kane, G.L. (ed.): Perspectives on super- 
symmetry II 1-153 (1997). 
[9] C. Gross, G. Marques Tavares, C. Spethmann, and 
M. Schmaltz, arXiv: 1209.6375 [hep-ph] (2012). 

[10] G. Aad et al. (ATLAS Collaboration), Eur. Phys. J. C 
71, 1828 (2011). 

[11] S. Chatrchyan et al. (CMS Collaboration), 
arXiv: 1302.0531 [hep-ex] (2012). 

[12] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 
71, 032001 (2005). 

[13] C. S. Hill, NucL lustrum. Methods A 530, 1 (2004). 

[14] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 
45, 001448 (1992). 

[15] CDF uses a cylindrical coordinate system with the z axis 
along the proton beam axis. For a particle or a jet, pseu- 
dorapidity is 77 = — ln[tan(0/2)], where 6 is the polar an- 
gle relative to the proton beam direction, and (j> is the az- 
imuthal angle while transverse momentum pt ~ \p\ sinO, 
and the transverse energy Et ~ E sin 9. 

[16] A. Bhatti et al, Nucl. Instrum. Methods A 566, 375 
(2006). 

[17] J. Alwall, P. Demin, S. de Visscher, R. Frederix, M. Her- 

quet, F. Maltoni, T. Plehn, D. L. Rainwater, and T. 

Stelzer, J. High Energy Phys. 09 (2007) 028. 
[18] J. Pumplin et al. (CTEQ Collaboration), J. High. Energy 

Phys. 07 (2002) 012. 
[19] T. Sjostrand et al., Comput. Phys. Commun. 238, 135 

(2001), version 6.422. 
[20] E. Gerchtein and M. Paulini, |arXiv:physics/0306031| 

(2003). 

[21] M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau, and 

A. D. Polosa, J. High Energy Phys. 07 (2003) 001. 
[22] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 

84, 052010 (2011). 
[23] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. D 

73, 032003 (2006). 
[24] A. Read, Nucl. Instrum. Methods A 425, 357 (1999). 
[25] A. Read, J. Phys. G: Nucl. Part. Phys. 28, 2693 (2002); 

T. Junk, Nucl. Instrum. Methods A 434, 425 (1999). 
[26] A. D. Martin, W. J. Stirhng, R. S. Thorne, and G. Watt, 

Eur. Phys. J. C 63, 189 (2009). 
[27] W. Beenakker, M. Kramer, T. Plehn, M. Spira, and 

P. M. Zerwas, Nucl. Phys. B515, 3 (1998). 



